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B 

Q , De Haas-van Alplien oscillations are studied for Fermi surfaces illustrating 

CJ I the Pippard's model, commonly observed in multiband organic metals. Field- 

and temperature-dependent amplitude of the various Fourier components, 
. linked to frequency combinations arising from magnetic breakdown between 

(T^ . different bands, are considered. Emphasis is put on the Onsager phase fac- 

J^ ! tor of these components. It is demonstrated that, in addition to the usual 

\Q I Maslov index, field-dependent phase factors must be considered to precisely 

— u ■ account for the data at high magnetic field. We present compelling evi- 

O ' dence of the existence of such contributions for the organic metal ^-(BEDT- 

TTF)4CoBr4(C6H4a2). 
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1. Introduction 

Fermi surface (FS) of numerous organic metals is an illustration of the 
textbook model proposed by Pippard more than fifty years ago to compute 
Landau band structure induced by magnetic breakdown (MB) in multiband 
metals [l|. This is the case of the FS of the strongly two-dimensional 
charge transfer salt ^-(BEDT-TTF)4CoBr4(C6H4Cl2) (where BEDT-TTF 
stands for the bis-ethylenedithio-tetrathiafulvalene molecule), which is re- 
ported in Fig. [1] [2| . Organic metals with such a FS are known to give rise to 
magnetic oscillations spectra involving linear combinations of the frequencies 
linked to the basic orbit a and the MB orbit /9. These frequencies correspond 
not only to MB orbits such as /9 -|- a or harmonics but also to 'forbidden fre- 
quencies' such as /3 — a that are not predicted by the semiclassical model of 
Falicov-Stachowiak [3|, Ij]. 

Only recently, analytic tools have been provided to account for the field 
and temperature dependence of the Fourier amplitude relevant to the various 
frequencies observed [2] • Besides, to our knowledge, little attention has been 
paid to the Onsager phase factor, yet. Though, according to Slutskin and 
Kadigrobov |5| and Kochkin [6|], a field-dependent Onsager phase should be 
observed for the considered FS, provided the magnetic field is large enough 
compared to the MB field. Almost ten years later, the same result was 
independently derived |8| in order to account for the discrepancy between 
calculations, which are valid for the low field range, and the experimental 
data for the lens orbit of Cd [7| which share similarities with the a orbit 
of Fig. [H However, still to our knowledge, no further study in this field 
has been reported up to now. In order to address this question, this paper 
is focused on the Onsager phase factor of the various Fourier components 
observed in the de Haas- van Alphen (dHvA) spectrum of the organic metal 
^-(BEDT-TTF)4CoBr4(C6H4Cl2) in fields of up to 55 T. 

2. Model 

In this section, we first present the model accounting for the field and 
temperature dependence of the amplitude of the various Fourier components 
entering the oscillation spectra (2[. In the second step, the field-dependent 
Onsager phase is considered. 




Figure 1: (Color online) Fermi surface of 6l-(BEDT-TTF)4CoBr4(C6H4Cl2) in the 
extended zone scheme. Blue solid lines depict the first Brillouin zone. Pink lines display 
the classical orbits considered for the data analysis, and arrows indicate the quasi-particles 
path on the principal orbits a and l3. Blue circles and diamonds indicate the turning points 
in the direction parallel to the chains and the Bragg reflection points, respectively. 



2.1. Fourier amplitude 

As displayed in Fig. [H the FS is composed of the a quasi-two-dimensional 
closed tube and a pair of quasi-one- dimensional sheets separated from the a 
orbit by a gap liable to be overcome by MB. Numerous classical MB orbits 
can be defined {rj = a, (3, a + (3, 2/3 — a, 2/3, etc.), the area of which are 
linear combinations of those relevant to the a and /3 orbits. The area of the 
latter is equal to that of the first Brillouin zone (FBZ). It can be remarked 
that 2/3 corresponds to both the classical orbit displayed in Fig. [T]and to the 
2nd j^armonic of /3. 

To account for this FS, a two-band system with band extrema Ao(i) and 
effective masses t?2o(i) (in units of the electron mass me) is considered |2|. 
The band gives rise to the quasi-one-dimensional part of the FS of Fig. [1] 
whereas the a orbit is built on the band 1. Assuming parabolic dispersion, 
the relevant frequency is Fa = mi{fi — Ai). The /3 orbit, generated by four 
tunnelings at the junction points, is built on both bands and 1 and, still for 
a parabolic band, has a frequency corresponding to the first Brillouin zone 
area, F^ = ma{fi — A^) + mo(/i — Aq) = rriis^fi — A^). In this case we identify 
the mass mo + rria with the mass m^ of the orbit /3. 

To compute the oscillating part of the magnetization at fixed number N 
of quasi-particles, we need to consider the oscillatory part of the free energy, 
defined by 

F{T,N,B) = Q{T,fi,B) + Nfi (1) 

For a constant N, the oscillatory part of the grand potential Q for a 
sample slab with area A can be written 

+ ^Y.11 ;^^^^'^(^) cos{2npFjB + p^,)- (2) 

p>l ri ^ ^ 

Damping factors can be expressed as Rr,^p{B, T) = R^p{B, T)R^p{B)R^^p{B)W^p 
where: 



Rl^ = pX^siYAi-\pXr,), (3) 
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R^^^ = eM-puom^TDB-^), (4) 

K,,^ = {^VoTHq^r^ (5) 

R'^p = cos{iTg^m^/2). (6) 

The field-and temperature-dependent variable (X^) and the constant (mq) 
are expressed as X^^ = uquIj^T / B and Uq = 27i'^kBme{eh)^^ = 14.694 T/K. 
The tunnehng (p^^ and reflection (go) probabilities are given by Pq = e'-^^/ss 
and pI + Qq = 1 [4] . 0o = h/e is the magnetic flux quantum, T^ is the Dingle 
temperature defined by Td = h{27ikBT)~^, where r~^ is the scattering rate, 
Bq is the MB field, rrirj and Qr, are the effective masses and effective Lande 
factor, respectively. In the case where the magnetic field direction is not 
parallel to the normal to the conducting plane (angle 6), B is changed to 
Bcos6 and the spin damping factor is written RL = cos{ti g^nirj/ 2 co?, 9). 
For convenience, energies {E) such as /i, Aq, Ai are expressed in units of 
Tesla, using the conversion {uQ/kB)\T / J] x E[J] = E[T]. B and T are 
the magnetic field [T] and temperature [K], respectively. Effective masses 
are expressed in units of the electron mass rrie, and magnetization in Tesla 
units. The advantage of taking this convention is that field, frequencies and 
temperature are not expressed in reduced units. 

Frequencies Fr^[T] can be written as F^ = m^(yU — A^) and are dependent 
on the chemical potential fi since they are proportional to the area enclosed 
by the orbits. Coefficients C^ are the symmetry factors of orbits r]. Namely, 
Ca = Ci3 = C2i3-a = 1 and Ca+i3 = 6*2/3 = 2. Integers n^ and n^ are the 
number of MB-induced tunnelings and reflections, respectively, yj^ is the 
Onsager phase factor of the orbit 77, defined by the number of turning points, 
i.e. 'n 12 times the number of extrema of the orbit along one direction (see 
Fig. [1]). A^ is given by dVi/dii = —N, and the chemical potential satisfies the 
following implicit equation: 

D 1 Tp 

/^ = /^o y^Yl —^v^vAT) sm{2TTp-^ + pep,,), 

which can be rewritten as 



/^ = /io-— ^M,(S). (7) 



where /zq is the zero-field Fermi energy. For a compensated system, in 
which case A^ = 0, it is equal to /iq = ("^o'^o + i^a^a) / ij^a + i^a)- The 
oscillatory part of the magnetization is defined as 



"^°-[^] = -MS dB • (') 

In this expression, the free energy, after factorization and simplification, 
is given by 

+ f E E ^ -^ [<^ - 2-p5 E ^V(B) + P^,) + est 

Since oscillating factors M^ entering Eq. [7] are at first order in damping 
factors Rr,,p{B,T) and small compared to /xq, Eq. |H]is solved at second order 
in Rr^p{B,T) (the first order part corresponding to the Lifshitz-Kosevich 
(LK) semi-classical result). This leads, after some algebra, to an expansion 
in power terms of the amplitudes 



V P> 



= -EE5^^-(B.r)sm(2.p§+pJ 

r, p>l ^ '' V ^ / 

+ E E ^^^i^Kp{B.T)R,,AB.T) [sin('2vr^^^^+py,,+pV.') 
- sml27r '— + pf^r^ - p '^q' \ 



r],r]' p,p 

+ ••■ (10) 



where the next terms are third order. From this step onwards, frequencies 
Frf are evaluated at /x = /xq: -^jj = ''^7j(/^o ~ ^i?)- According to the above 
expression, magnetization spectrum can now be expressed in terms of both 
classical and non-classical frequencies, still noted as F^ in the following, and 
can be expanded as: 



rrir. 



Y^ Apr^ sin ( 27rp-^ + P(prj] ■ 



(11) 



It is important to stress that the amphtude Apr^ involves not only the 
contribution of the p*'^ harmonics of the rj classical orbit, given by the LK 
formalism (Ap^ oc Rr],p) but also higher order corrections, calculated here 
at the second order in damping factors. The expressions of the dominant 
Fourier components, considered for the data analysis, are given belowo 
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(13) 
(14) 



+ 2i?a,l-R2/3-Q:,l — -R/3,2-R2/3,2 " Ri3,aR2I3,1 
Fp-a 



Rji,! ~ q-R/3,l-R/3,3 ~ -jRl3,2RfS,4 ~ Ra,2Ra+fS,2 



irrriis 



[-Ra,l-R/3,1 + Ra,2Ra+^,l + Rl3,2Ra+l3,l + Rl3,lR2(^-a, 
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F2l3-a r, F2l3-a 



Ra,lRl3,l — '^Ra+P,2Ra+P,l — --R^,3-R2/3-«,l 



-R 



11171213-0 

F2l3-2a 

7117113 



2/3-a,l 



71-112(3 



1 1 

— Ra,lR/3,2 + -Ra,1,Ra+l3,2 



2-^0,2-^2/3,1 + 2-R2;9-a,2-R2^,l + 2-Ra, 1-^2/3-0,1 + -Ra,4Ra+l3,2 



1 1 

+ --Ra,2-R/3,2 + --R/3, 2-^2/3-0,2 

As discussed in Ref. |2i], the leading term of Eqs. [121 [HI and [HI relevant 
to the classical orbits, a, (3 and 2/3 — a, respectively, corresponds to the LK 



(15) 
(16) 

(17) 

(18) 



(19) 



^In Ref. I2I], all the Onsager phases are arbitrarily set as 0. For this reason, amplitudes 
in Eqs. [T2| [T4l and [TSl have an opposite sign compared to the data in Ref. [2|, accounting 
for the TT dephasing reported in Tabled] 
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Figure 2: Phase w as a function of the argument Bo/{2ttB). For smah fields B compared 
to the magnetic breakdown field Bq, the phase vanishes while it goes to 7r/4 at large fields. 



formalism. This statement also holds for /3 + a, 2a and 2/3 even though 
the second order terms may have magnitude close to the first order term, 
being able to yield non-monotonous field and temperature dependence. In 
the specific case of Eq. [T51 the leading term involves the contributions of 
both the classical orbit 2/3 displayed in Fig. [T] and the second harmonics of /3 
which are accounted for by the damping factors -R2/3,i and -R/3,2, respectively. 
In contrast, there is no first order term entering Eqs. [T6l and [T9] relevant to 
P — a and its second harmonics, respectively, since these Fourier components 
correspond to 'forbidden frequencies'. 

2.2. Onsager phase factor 

Turn on now to the determination of the Onsager phase factors (/)^ entering 
Eq. [m Besides a phase factor (p^, deduced from the classical orbits' phases 
appearing in Eq. [21 an additional phase factor u is added to (p^j each time 
a quasiparticle is reflected at a MB junction. Indeed, according to Refs. 
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p, la |8[ , the matrix for the incoming and outgoing wave- function amphtudes 
at each junction point is given by 



M = ( ^°^"^^ '^l ) (20) 



with 



TT R 

uj{B) = h a:;log(x) — X — argr(za;), x = — . (21) 

After a reflection, the quasi-particle amphtude takes a factor goexp(— iw) 
and goexp(ico') for quasi-particle path orientation clockwise and counter- 
clockwise, respectively. As displayed in Fig. [21 uj goes to zero at low field. In 
contrast, it takes noticeable values as the field is larger than Bq, going to 7r/4 
at large field. For example, for Bq = 35 T, and -B = 55 T, that are relevant 
values for the compound studied in this paper, we obtain u = O.lGvr, which 
is not negligible, especially if the number of reflections events nj^ is large. 
According to Eq. |20l the Onsager phase factor is given by 

(j)r, = ^r,- ^"v^^B) (22) 

The ifr^ and rC values relevant to the Fourier components appearing in 
Eqs. [12] to [19] are given in Table [1] We notice in particular that the index 
rf^ can be negative, due to algebraic combinations of the individual phases 
present in the sine function of Eq. [10] It can also be remarked that the 
Fourier component with frequency F2/3 arises from the second harmonics of 
/3 and the 2/3 orbit displayed in Fig. [T] Nevertheless, these two contributions 
have the same Onsager phase. Besides, for a given 77 Fourier component, all 
the involved second order terms (see Eqs. [12] to [T9|) can be viewed as arising 
from algebraic combinations of classical orbits yielding the same Onsager 
phase. 

3. Experimental 

Field- and temperature-dependent magnetic torque of the considered crys- 
tal was studied in Ref. [2|. It was synthesized by electrocrystallization tech- 
nique as reported in Ref. (9|. Its size is approximately 0.12 x 0.1 x 0.04 mm^. 
Recall that magnetic torque was measured with a commercial piezoresistive 
microcantilever, in pulsed magnetic fields of up to 55 T with a pulse decay 



Table 1: Onsager phase factors predicted by Eq. [22]for the various Fourier components 
considered for the data analysis and (/3,,/7r values accounting for the best fits of Eqs. [12] 
to [19] to the oscillatory torque data in the temperature range 1.4 - 4.2 K (see Fig. [3] 
for the data at 1.4 K), assuming either constant (assumption (i): n^^ — in Eq. [22]) or 
field-dependent Onsager phase factor (assumption (ii)). 



Fourier 


Predicted 


(Eq.[22D 


fittings 


0r, = V^r, 


fittings 0r; = 


<^r, - n^W 


component 


ip^JTl 


n'n 


F(T) 


(/^^/vr 


F(T) 


V^^/vr 


a 


1 


2 


949.7 ± 1.2 


0.75 ± 0.04 


947.3 ± 1.4 


1.16 ±0.05 


2a 





4 




-0.52 ± 0.21 




0.30 ±0.13 


P-a 





-2 




0.35 ± 0.21 




-0.06 ± 0.21 


/3 


1 





4631 ± 6 


1.08 ± 0.19 


4631 ± 6 


1.08 ±0.19 


/3 + a 





2 




-0.7 ± 0.8 




-0.3 ± 0.8 


2/3 -2a 





-4 




0.5 ± 0.7 




-0.3 ± 0.7 


2(5 -a 


1 


-2 




-0.54 ± 0.34 




1.05 ± 0.34 


2/9 


1 







1.1 ± 0.5 




1.1 ±0.5 



duration of 0.32 s. Variations of tlie cantilever piezoresistance were measured 
in tlie temperature range from 1.4 K to 4.2 K witfi a Wlieatstone bridge witli 
an ac excitation at a frequency of 63 kHz. Tlie angle between the normal to 
the conducting plane and the magnetic field direction was 9 = 1°. 

4. Results and discussion 

The field-dependent amplitudes A^ of the various Fourier components 
entering the spectra can be analyzed on the basis of Eqs. [12] to [19], keep- 
ing in mind that they are related to torque oscillations amplitudes AL as 
A^ = roA^/(i?tan^) where tq is a prefactor depending of the cantilever stiff- 
ness, crystal mass, etc.. Onsager phase factors are considered within the 
assumptions of (i) constant (n^ = 0, i.e. 0^ = ipr^ in Eq. [22]) and (ii) field- 
dependent 0^. 

The field range above 20 T, in which dHvA oscillations are observed is 
considered for the data analysis. As discussed above, numerous physical 
parameters enter the oscillatory spectra. In order to reduce the number of 
free parameters for the fittings, effective masses [rria = 1.81, m/j = 3.52), MB 
field {Bo = 35 T) and Dingle temperature (To = 0.79 K for all the classical 
orbits) are taken from Ref. |2|] and kept fixed. Besides, the effective Lande 
factor is assumed to be the same for all the orbits ((7^ = g*). As a result, 
besides prefactors Tq and the Lande factor g*, only the two frequencies F^ 
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Figure 3: (Color online) Oscillatory torque data at 1.4 K (solid circles) and best fit of 
Eqs. [12] to [19] (green solid line). The residuals (open squares) and the various Fourier 
components entering the fittings (solid lines) are shifted down from each others. 
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and Fp and the various field-independent parts of the Onsager phase factors 
ifrj are free parameters. 

As an example, the high field range of the data measured at 1.4 K are 
reported in Fig. [31 Best fits obtained within either assumptions (i) or (ii) 
are indiscernible to the naked eye in the field range considered for the data 
analysis, even though, strictly speaking, oscillations are no more periodic in 
1/B within Eq. [22J Actually, it can be checked that, due to the limited field 
range in which oscillations are observed, Fourier transforms of the best fits 
obtained within either the assumption of (i) constant or (ii) field-dependent 
Onsager phase, are indiscernible, as well. The deduced effective Lande factor 
is g* = 1.9 ± 0.2, in agreement with the reported value of Ref. (2|. Deduced 
values of ^9^, F^ and F^ are given in Tabled! 

As expected, Fp, ipj^ and ip2p are insensitive to the considered assumption 
(i) or (ii) since only tunnelings enter /3 and 2/3 orbits {np = and n2/3 = 
in Eq. |22|) . In addition, ipp is in agreement with the predicted value within 
the error bars. This is also the case of ip2i3 although a large uncertainty is 
obtained due to the small amplitude of this Fourier component (see Fig. [3]). 

In contrast, the value of both the frequency F^ and Onsager phases of the 
Fourier components involving a depend on the considered assumption (i) or 
(ii). Crudely speaking, the Fa value deduced from fittings within assump- 
tion (ii) accommodates to compensate the field-dependent phase. However, 
the observed change is small and remains within the uncertainty. Not any 
of these Onsager phase values, deduced within assumption (i), are in agree- 
ment with the predictions (see Tabled]). In other words, field- independent 
Onsager phases given by it/2 times the number of orbit extrema in fc-space 
cannot account for the data. Assuming a field-dependent phase factor, ipa 
and ip2a come much closer to the predictions of the model, even though the 
discrepancy with the predicted values are still slightly off the error bars. Be- 
sides, ipi3-a and ip2i3-a are in very good agreement with the predictions. This 
statement also stands for v^/^+a and ^2p-2a albeit the error bars are very large 
due to the small amplitude of these Fourier components. 

5. Conclusion 

Many years ago, theoretical calculations predicted that Onsager phase 
factor of quantum oscillations includes a field- dependent part in the case of 
a magnetic breakdown orbit with refiections J5|, |6|, |8|. To our best knowl- 
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edge, this feature was only considered in the case of Cd [7| which is a three- 
dimensional elemental metal. 

It is demonstrated that field-dependent phase is necessary to account 
for the oscillatory spectrum of the two-dimensional organic metal 6'-(BEDT- 
TTF)4CoBr4(C6H4Cl2). The Fermi surface of this compound achieves the 
network of coupled orbits model proposed by Pippard more than fifty years 
ago[l| which is relevant for many organic compounds. In agreement with 
the above mentioned theoretical predictions, magnetic breakdown orbits in- 
volving reflections, namely all the orbits including the a component for the 
considered Fermi surface topology, are accounted for by field-dependent On- 
sager phase factors. This result confirms that field-dependent phase factor is 
a general feature of magnetic breakdown orbits. 

As a result, the value of Fq, deduced within Eq. [22] is slightly reduced com- 
pared to the value derived assuming constant Onsager phases (i.e. through 
either direct fitting assuming 0^ = y9^ or Fourier analysis of the data), even 
though the observed discrepancy stays within the error bars. Despite field- 
dependent phase factors, the magnetic oscillations periodicity in XjB is pre- 
served, within the experimental uncertainty. This behaviour is likely due to 
the relatively small field range considered (20 T to 55 T). Obviously, keeping 
in mind that the magnetic breakdown field of the studied compound is rather 
large (B^ = 35 T), larger effects are expected for compounds with smaller 
magnetic breakdown field. 
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